home *** CD-ROM | disk | FTP | other *** search
/ CD ROM Paradise Collection 4 / CD ROM Paradise Collection 4 1995 Nov.iso / os2 / adaptor.zip / ADAPT.ZIP / adaptor / examples / f77 / readme < prev   
Text File  |  1993-03-23  |  1KB  |  51 lines
  1. Programs contained in this directory contain no array operations (only
  2. using FORALL or parallel loops).
  3.  
  4. Therefore the flag -F77 or -F90 have the same meaning (F77 loops will
  5. not translated to F90 array operations).
  6.  
  7. Both programs check their results by comparing the results with the 
  8. original values.
  9.  
  10. a) FFT
  11. ======
  12.  
  13. fft.f: Fast Fourier Transformation (indirect addressing)
  14.  
  15.   Algorithm applies to a random sequence the fft and its inverse. It
  16.   compares the original sequence with the new computed sequence.
  17.  
  18.   Output (should be near to zero):
  19.  
  20.       Max =   2.2942879E-07
  21.  
  22.   On Alliant FX/2800 there are problems with optimizations switched on
  23.   (use -Og instead of -Oguv)
  24.  
  25.   With Adaptor 1.0 the old version runs in some problems due to the array
  26.   normalization (lower bounds of arrays must be 1).
  27.  
  28. b) LU and MATMUL
  29. ================
  30.  
  31. lumm.f: LU Decomposition and Matrix Multiplication
  32.  
  33.   The Algorithm takes a n x n matrix with random values and makes an
  34.   LU Decomposition. Then it computes the inverse to a.
  35.  
  36.   Finally a and its inverse are multiplied. The result must be the
  37.   unity matrix.
  38.  
  39.   Output (maximal eps should be near to zero):
  40.  
  41.      Arrays are allocated
  42.      PROGRAM STARTS
  43.      LU ready
  44.      now finds Y with L * Y = E
  45.      now finds X with R * X = E
  46.      b is the inverse
  47.      lu = a * b
  48.      matmul correct, now check = E
  49.      maximal eps =   2.4130940E-04
  50.  
  51.